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Abstract 


Ay matrix element of four-quark operator relevant to its lifetime is rean- 
alyzed by QCD sum rule. The new ingredients introduced are that (1) 
more vacuum condensates are considered; (2) different quark-hadron dual- 
ity is adopted; and (3) the possible deviation from the vacuum saturation 
assumption for the four quark condensates is considered. With ку = 4 the 


related hadronic parameter т is calculated to be (3.6 + 0.9), and the lifetime 
ratio T(A,)/7(B°) = (0.83 + 0.04), which is consistent with the experimental 
data. 
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I. INTRODUCTION 


Heavy baryon lifetimes provide testing ground for the standard model, especially for QCD 
in some aspects, because they can be calculated systematically by heavy quark expansion 
[I]. Theoretically, if we do not assume the failure of the local duality assumption, the heavy 
hadron lifetime differences appear, at most, at the order of 1 / mi B]. Recent experimental 
result on the lifetime ratio of the A, baryon and B meson is [| 


т(Ль) 


which has shown some deviation from the theoretical expectation. This is one of issues 
concerned with the tests of the heavy quark expansions and has drawn a lot of theoretical 
attention HHB], and may imply the potential importance of the O(1/m;j) effect for the above 
heavy baryon and heavy meson lifetime difference ДӨ]. 

It is at the order 1/mg that appear four-quark operators, whose contributions to decay 
widths are enhanced by a phase-space factor of O(1677) with respect to the leading terms 
in the operator product expansion. Consequently one should include them in the prediction 
for non-leptonic decay rates [Ц], To the order of 1/mj, the lifetime ratio was calculated as 
follows |l]; 


т(М) Êl) — (B) eB) 
wm)" P ag cde 


tipi Bi (my) + po Bo(my) + рзє (ть) + paeo (ma) + [ps + ps B(my)]r(my)] » (2) 


where the term proportional to & = (f5/200MeV)? arises from the contributions of four- 
quark operators. The numerical values of the coefficients сс and p;'s (i = 1 — 6) have 
been calculated in Ref. Д. u2(H,) and џ2(В) are the averages of b-quark kinetic and 
chromomagnetic energy, ТӨ To the order of 1/m?, the lifetime ratio is 0.98. At 
the scale my, the values of p;'s are pı = —0.003, po = 0.004, рз = —0.173, p4 = —0.195, 
ps = —0.012, pg = —0.021. В, Bo, є, ө, r and B are the parameterization of the hadronic 
matrix elements of the following four-quark operators, 


(Bl, (1 — 75)997"(1 — 95)0|B) = By fam , 
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In Eq. (B), the energy scale for the parameters is at ть. In Eq. (B), generally the renormal- 
ization scale is arbitrary, and the parameters depend on it. It can be taken naturally at the 
low hadronic scale to apply the heavy quark expansion. 

These parameters should be calculated by some nonperturbative QCD method. The 
QCD sum rule [J], which is regarded as a nonperturbative method rooted in QCD itself, 
has been used successfully to calculate the properties of various hadrons. In Refs. [ and 
0]. the mesonic parameters B; and є; have been calculated by the QCD sum rules within 
the framework of heavy quark effective theory (HQET) [13]. The baryonic parameters r 
and B have been calculated in Ref. []. As a result, to the order of 1/ mj, the theoretical 
calculation still cannot explain the experiment result {). 


II. THEORETICAL CONSIDERATION 


In this paper, we try to match the experimental data by considering more subtleties in 
the theoretical analysis. As analyzed in Ref. [14], the lifetime ratio ([l] depends crucially 
on the value of the baryonic parameter r. Therefore, we shall carry out a reanalysis for 
baryonic parameter r by the QCD sum rule. We note that in the analysis of Ref. [|], more 
condensates could be included in. They are gluon condensate and six-quark condensate. It 
may be important because the dimension of the gluon condensate is four. In spite of being 
dimension nine, the six-quark condensate may also have significant contribution, because 
the diagram is not suppressed by any loop factor. ТИЯ. The from the heavy baryon masses 
tells us that the pure tree diagram is important [1 . Therefore in our calculation, these 
two kinds of condensates will be considered, in bd to those considered in Ref. [7]. 

The conclusion of B = 1 does not change. It follows from the valence quark approxima- 
tion which is used both in Ref. [ and in our QCD sum rule analysis. 

To calculate r, the following three-point Green's function is constructed, 


по) =? f dedye 7" (op; (r0 (0) 0010) . (5) 
where w = v - k and w! =v-k’. The Ag baryonic current j" was given in Refs. 19-9, 


J? = Cla +b orh , (6) 


where a and b are certain constants which will be discussed later, h, is the heavy quark field 
in the HQET with velocity v, C is the charge conjugate matrix, т is the flavor matrix for 


Ло, 
1 0 1 


In Eq. (8), Ó denotes the four-quark operator 


me c ns 
Ó- bw ; —P м (8) 


Note < A,|O|A, >= — < A,|O|A, > in the valence quark approximation F], where 
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In terms of the hadronic expression, the parameter r appears in the ground state contribution 
of (w, ш”), 


fMAglOIAg) 1+ 4 
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+ higher states , (10) 


where the " higher states" denotes the contribution of resonances and continuum. (Ag|O|Ag) 


has appeared in Eq. (H), A = ma, — то and the quantity fa is defined as 


(0"|Ag) = fau , (11) 


with u being the unit spinor in the HQET. The QCD sum rule calculations for fa were 
given in Refs. [T5] [IS]. On the other hand, this Green's function II(w,w') can be calculated 
in terms of quark and gluon language with vacuum condensates. This establishes the sum 
rule. For the resonance part of Eq. (10), we adopt the assumption of quark-hadron duality. 
The calculation of П(о, cw") in HQET is straightforward. The fixed point gauge [19] is 
used. The tadpole diagrams in which the light quark lines from the four-quark vertex are 
contracted have been subtracted. The following values of the condensates are used, 


(qq) = —(0.23 GeV)? , 
(a,GG) = (0.075 + 0.015) беу“, 
(990 wG"" gq) = то(а0) , то = 0.8 GeV? . (12) 


Note we have adopted the new value for the gluon condensate from the recent analysis of 
the heavy quarkonium spectrum [BO]. However, our final results are not very sensitive to 
the gluon condensate since its contribution is only about 6% of the whole sum rule as shall 
see later. 

Except for the quark gluon mixed condensate, our calculation would be consistent with 
that of Ref. [fj] if the gluon and six-quark condensates were omitted. While the calculation 
can be justified if (—w) and (—w’) are large, however the hadron ground state property should 
be obtained at small (—w) and (—w’). These contradictory requirements are achieved by 
introducing double Borel transformation for w and ш”. It is defined as 


B- lim lim cere | d y= (s) | (13) 
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There are two Borel parameters 7 and 7’. They appear symmetrically, so 7 = 7' = 2T is 
taken in the following analysis. The reason for the factor 2 is similar to that explained in 


the Ref. B1]. 


III. DUALITY ASSUMPTION 


Generally the duality is to simulate the higher states by the whole quark and gluon 
contribution above some threshold energy we. The whole contribution of the three-point 
correlator II(w, о”) can be expressed by the dispersion relation, 


Dou) = [^ av [^ cav) —— | (14) 


(v! — ш”) 


With the redefinition of the integral variables 


„+! 
V === 
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the integration becomes 


|" dv |" ПВ а = Ar dv. ре ы. (16) 
0 0 0 iy, 


It is in v4 that the quark-hadron duality is assumed [13], 


2 үә v ImII(v, и 
higher states — — dv, t jp UA . 
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(17) 


This kind of assumption was suggested in calculating the Isgur- Wise function in Ref. and 
was argued for in Ref. B2]. As pointed out in [P2[I3], in calculating three-point functions 
the duality is valid after integrating the spectral density over the "off-diagonal" variable 
и = i(v-v). 

The sum гше for (Ag|O|Ag) after the integration with the variable v_ is obtained as 


b A b? b 
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Ki (qa) V^) — ka T (18) 
where ку, K2 are the scale parameters used to indicate the possible deviation from the 
factorization assumption for the four- and six-quark condensates. куо = 1 corresponds to 
the vacuum saturation approximation. ку = (3 ~ 8) is often introduced in order to include 
the nonfactorizable contribution and to fit the data [B3]. There is no discussion of кә in 
literature so we use кә = 1. 

As argued in [ГЇ], the choice a = b = 1 tends to yield the optimal interpolating current 
for the A, baryon in HQET. We shall adopt a — b — 1 in our numerical analysis. 

The parameters f4, Л etc were obtained by the HQET sum rule analysis of two-point 


correlator [518], 


(a+b)? ,. 4 we jp a? + BP 2ab(qq) me 
SL fhe? = n бети? — TM (а _ Ty. (19) 


We have not included a, corrections in Eq. (9), because they are also neglected in the 
sum rule for r (189). The values of the parameters are f? = (2.9 + 0.5) x 107? GeV?, 
A = (0.9 + 0.1)GeV with the threshold we to be (1.2 + 0.2)GeV and the Borel parameter Т 
in the window (0.2 — 0.4) GeV. 

In order to minimize the dependence of the parameters we divide Eq. (18) by Eq. (19) to 
extract (Ag|O|Ag). The variation of the matrix element with we for two values of ку = 1,4 is 
given in Fig. 1. The value of w, is (1.2+0.1) GeV. The sum rule window is T = (0.15 —0.35) 
GeV, which is almost the same as that in the two-point correlator sum rule. We obtain for 
Кі = 4 


(Ag|O|Ag) = (1.6 + 0.4) x 107?GeV? . (20) 


By taking fp = 200 MeV, 


r = (3.6 + 0.9). (21) 

If we use ку = 1, we get 
(Ag|O|Ag) = (5.5 + 1.0) x 107?GeV? , (22) 
r = (1320). (23) 


In our numerical calculation, the contribution of the perturbative term is smaller than 
those of the condensates due to small T' in the duality window and high dimensionality of the 
spectral function ( v? in eq. (18) ). There are several known examples where the perturbative 
part has small contribution in the sum rule 1681104, without affecting the reliability of 
the results. Nevertheless, a hierarchical structure among various condensate terms exists in 
those cases and the present case. In the present case, because of the cancellation between the 
contributions of the quark condensate and the mixed condensate, as well as the smallness 
of the contributions of the gluon and the six-quark condensates, the influence of the four- 
quark condensate is significant and is enhanced by assuming the deviation from the vacuum 
saturation. 

In Ref. Д the authors used the following duality assumption, 


1 г со П 4 
higher states = — dv | dp DNA , (24) 
T Jwe we (v — w)(v' — ш) 


In their calculation the matrix element (Ag|O|Ag) increases from 0.4 x 10-?GeV? to 1.2 x 
107?GeV? in the working region of the Borel parameter when we varies from 1.1 GeV to 
1.3 GeV. In other words, their analysis is very sensitive to the continuum threshold, which 
might imply the above duality assumption is not good. In contrast, the dependence on w, 
is not so strong in our approach. 

The 1/mij corrections to the above results can be analyzed in principle. While having 
little influence on our above calculation, they formally belong to the O(1/m) effects to 


6 


the A; lifetime. It should be noted that the value of r we have obtained above is at some 
hadronic scale, other than the scale ть, because we have been working in the framework of 
the HQET, in which the natural scale is [nag K ть. The renormalization group evolution 
of the relevant operators was calculated in Ref. [25,4]. Information on parameters В andr 
at scale rm is necessary to obtain the lifetime ratio of Eq. (3. By choosing &s(Hhaa) = 0.5 
(corresponding to Hhaa ~ 0.67 GeV), Ref. HJ gives 


т(ть) = [1.54 + 0.18B (haa) |? (lhaa) н 


В(ть) = tte f (25) 
That is 
В(ть) ~ 0.58, (26) 
т(тъ) = (8.2 + 1.6) (27) 
for ку = 4 or 
r(m,) = (2.3 + 0.6) (28) 
for ку = 1 from Eqs. (21) and (63). 
The A, and B? lifetime ratio given in Eq. (Q) is expressed specifically as 
A = 0.98 — 0.17e; (my) + 0.20е›(ть) — (0.013 + 0.022 B (mj))r(m;) 
= (0.83 + 0.04) (29) 
for ку = 4 or 
BE = (0.93 + 0.02) (30) 
for kı = 1. Where the values єү(ть) = —0.08 and €2(m,) = —0.01 have been taken from 


the QCD sum rules [E]. From Eq. (BO), we see that with the vacuum saturation (кү = 1), 
although r is enhanced by about six times compared to that in Ref. ff, it is still not large 
enough to account for the data Eq. (1). The life time ratio between A, and B mesons can 
be explained if we also take into account the nonfactorizable contribution of the four-quark 
condensate, as can be seen from Eq. (29). 


IV. SUMMARY AND DISCUSSION 


In summary, we have reanalyzed the QCD sum rule for the A, matrix element of the four- 
quark operator relevant to the lifetime of 4j. Compared to the previous analysis [/], the new 
ingredients we have introduced are that (i) more condensates are considered; (ii) different 
quark-hadron duality is adopted and (iii) the possible deviation from the vacuum saturation 


T 


assumption for the four quark condensates is considered. Of these ingredients, the last two 
points have more important influence on numerical results than the first point by noting 
the effect due to more condensates is small (~ 10%). The second point is more essential 
than the third point. Note that due to the second point, r is about six times enhanced by 
comparing the result of Eq. (23) with that in Ref. [7]. And due to the third point, r is about 
three times enhanced by comparing the results of Eqs. (23) and (61). With fixed duality 
assumption for the sum rule, our result shows that the baryonic parameter r is significantly 
dependent on the nonfactorizable effect of the four-quark condensate and r(m;) = (2.3+0.6) 
for kı = 1 and т(ть) = (6.2 1.6) for ку = 4. With this latter value the A, baryon and B? 
meson lifetime ratio has been calculated to be 7(A,)/7(B°) = (0.83 + 0.04), which is close 
to the experimental data 7(A,)/7(B°) = (0.79 + 0.06). From the above results we draw the 
conclusion that under the local duality assumption it is possible to resolve the issue of the 
A, baryon and B? meson lifetime ratio in the framework of heavy quark expansion. 
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Figure captions 
Fig. 1. Sum rule for (Ag|O|Ag), where O has been given in Eq. (8). From top to 


bottom these curves correspond to we = 1.3, 1.2, 1.1 GeV respectively. The sum rule window 
is T = (0.15 — 0.35) GeV. 
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